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1  Accomplishments 

The  various  lines  of  progress  were: 


1.1  Nonlinear  Plants 

Over  the  last  decade,  the  extension  of  H°°  control  to  a  nonlinear  system  has  been  converted 
to  solving  two  particular  PDE’s.  One  partial  differential  inequality,  P.D.  I.,  the  state  feedback 
HJBL  corresponds  to  the  CT RL 00  problem  where  the  controller  can  access  the  full  state  of 
the  plant.  Indeed  solving  this  PDE  gives  an  optimal  solution  to  the  state  feedback  CTRL 00 
problem  and  can  be  solved  off  line.  The  second  PDE,  the  information  state  PDE ,  gives  the 
dynamics  of  the  controller  and  so  it  must  be  solved  ON  LINE. 

Unfortunately,  these  are  P.D.E’s  on  the  state  space  of  the  original  plant  (often  a  high  di¬ 
mensional  space)  so  a  numerical  solution  faces  what  is  called  the  curse  of  dimensionality. 
Getting  around  this  is  the  main  challenge  to  the  field  and  it  requires  basic  theory  as  well  as 
numerical  analysis. 

1.1.1  State  Estimation 

Collaborators  and  I  have  done  considerable  work  on  both  of  these  equations.  Most  of  the  ef¬ 
fort  goes  to  the  information  state  PDE,  since  at  first  glance,  computation  (on  line)  seems 
impossible  in  general  situations  with  state-space  higher  than  two.  However,  recent  mathe¬ 
matical  work  with  M.  James  and  partly  with  Bill  McEneaney  and  Peter  Dower  discovers  that 
in  various  important  cases  the  controller  dynamics  in  practice  might  NOT  suffer  prohibitively 
from  the  curse  of  dimensionality.  One  point  is  that  the  information  state  is  often  supported 
on  a  0,1,  or  2  dimensional  manifold  M  and  then  the  computation  required  to  propagate  it  is 
manageable.  Thus  pinpointing  such  cases  is  important  though  mathematically  challenging. 
Cases  under  current  investigation  are 

(a)  (with  James)  [HJ99book,prep] 

Mixed  sensitivity=  2-block,  with  dim  M  is  the  antistable  manifold  of  the  plant. 

(b)  (with  James,  McEneaney  )  [cdc98] 

Measure  all  but  k  states  directly  (cheap  sensor  case)  Example: 
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1.  Missile  with  no  air  probes  k  =  number  of  missing  airprobes.  dim  M  =  k 

2.  Compressor  dim  M.  =  1 

If  #  “noise  free  sensors”  <  #  disturbances,  then  dim  M  =  0. 

This  work  has  spread  over  several  years  with  the  main  additions  from  the  last 
year  giving  reassuring  estimates  demonstrating  how  small  amounts  of  noise  do  not 
radically  change  results. 

(c)  (with  Dower,  James)  [DHJacc99]  Linear  systems  with  non-linear  actuators-  We  examine 
an  obvious  architecture  for  controllers:  just  feed  the  actuated  control  signal  into 
the  state  estimator  (rather  than  the  control  signal).  This  removes  the  direct  effect 
called  ’’windup”.  If  the  noise  level  is  low,  then  the  information  state  equation 
reduces  to  an  ODE,  so  it  is  solvable.  This  is  a  dim  M  =  0  situation.  Even  if  there 
is  substantial  noise,  our  theory  of  nonlinear  control  is  broad  enough  to  prove  in 
many  cases  that  this  construction  is  the  best  that  can  be  done. 

The  practical  interpretation  of  the  constructions  above  in  (a)  and  (b)  are  that  they  tell  us 
where  to  put  grid  points  in  numerical  computation.  Ideally  all  grid  points  go  on  the  (small) 
manifold  M..  This  is  risky.  In  practice  one  should  place  some  points  near  M.  and  a  few 
far  from  A4.  The  appropriate  mathematics  for  determining  placement  of  grid  points  off  of 
M  appears  to  be  s  perturbations  of  solutions  to  the  information  state  equation.  This  falls 
into  the  framework  of  singular  perturbation  theory,  and  is  the  next  generation  of  work  along 
these  lines. 


1.1.2  State  Feedback 

Our  work  on  the  state  feedback  HJB  equations  focuses  on  solving  them  numerically. 
Work  is  with  Michael  Hardt-Ken  Kreutz  Delgado,  ECE  Dept.  These  are  old  equations  and 
our  approach  is  to  take  the  most  refined  software  publically  available  (which  originally  was 
developed  for  path  planning)  and  convert  it  for  our  purposes.  There  was  gratifying  success 
with  the  biped  walking  problem  below.  To  set  it  in  context,  we  mention  the  compressor 
problem  completed  within  the  last  year  though  focused  on  the  year  before. 

(a)  Moore  Greitzer  Compressor  Stall  Model:  This  has  a  3  dimensional  statespace  with 

saturation  constraints,  and  our  technique  seems  to  work  fine.  Ultimately,  it 
produces  an  explicit  26  parameter  (not  so  big)  neural  net  control  law.  The 
technique  should  extend  to  5  dimensions  with  little  modification. 

(b)  Classical  Biped  Walking  Problem-  full  5  link  model:  This  is  a  horrendous  path  plan¬ 

ning  as  well  as  modeling  problem.  The  state  space  is  14  dimensional  with  two 
“phases”.  One  foot  on  ground  imposes  2-algebraic  constraints,  two  feet  on  the 
ground  imposes  4-algebraic  constraints.  There  is  saturation  of  actuators  and 
states.  Finding  an  energy  minimizing  motion  amounts  to  solving  a  HJB  at  one 
point  in  state  space.  Work  in  the  literature  is  with  vastly  compromised  models, 
so  we  think  our  solution  to  this  is  a  considerable  advance. 

Item  (b)  gives  further  evidence  of  the  soft  underbelly  of  the  Bellman  Equation.  Given  a 
single  point,  even  in  a  10  dimensional  space,  one  can  compute  the  value  function  at  that 


point.  (Note  this  is  in  gross  distinction  to  elliptic  equations).  The  main  difficulty  is  less  in 
solving  the  Bellman  equation,  than  it  is  in  determining  where  you  want  to  solve  it  and  in 
simplifying  data  on  an  irregular  grid  in  high  dimensional  space. 


1.2  Symbolic  Algebra  for  Systems  Research 

We  are  the  providers  to  Mathematica  of  general  Noncommuting  algebra  capability.  We  have 
a  major  computer  package,  NC Algebra,  which  is  a  collection  of  ’’functions”  for  Mathematica 
designed  to  facilitate  manipulation  and  reduction  of  noncommutative  algebraic  expressions 
in  operator  theory  and  engineering  systems.  Noncommutative  algebra  is  a  cornerstone  of 
matrix  theory  and,  consequently,  of  research  in  linear  engineering  systems.  It  also  occurs  in 
many  other  fields. 

Our  work  has  various  phases  which  run  from  bread  and  butter  to  ambitious  research.  NC  Al¬ 
gebra:  The  basic  commands  of  Mathematica,  Maple  etc  do  NOT  apply  to  noncommuting 
expressions.  Our  package  NCAlgebra  does  basics  of  noncommutative  algebra  in  Mathe¬ 
matica.  Simplifying  expressions.  Eliminating  unknowns  from  collections  of 
equations:  This,  the  focus  of  current  efforts,  is  the  strict  noncommuting  generalization 
of  the  commutative  Solve  and  Eliminate  commands  of  Mathematica  or  Maple.  It  uses  our 
NCGB  code  plus  algorithms  of  our  own  construction  for  eliminating  redundant  equations. 

A  major  focus  of  research  is  using  and  extending  these  algorithms  to  help  an  engineer  ’’dis¬ 
cover”  formulas  which  are  critical  to  his  particular  problem.  This  is  ambitious  but  the  stakes 
are  high.  We  do  many  experiments.  We  are  steadily  building  a  collection  of  classic  theo¬ 
rems  and  formulas  which  can  be  ’’discovered  semi-automatically”;  this  is  how  we  judge  our 
methods.  Also  we  make  new  discoveries. 

Our  work  indicates  that  our  software  is  useful  in  the  area  of  singular  perturbations.  There 
one  gets  truly  long  and  unpleasant  formulas  which  a  human  manipulates  slowly  and  with 
uncertain  accuracy.  Our  computer  methods  are  accurate  and  our  NCGB  plus  cleaning  plus 
sorting  algorithms  could  well  save  the  user  considerable  amounts  of  time.  In  classic  cases  our 
software  produces  the  classical  equations.  In  other  problems  being  studied  with  McEneaney 
(see  above)  we  find  this  symbolic  software  a  valuable  tool. 

My  student  D.  Kronewitter  just  determined  the  solution  structure  of  a  class  of  3  x  3  matrix 
completion  problems.  There  were  31,000  cases  modulo  about  500  permutations  which  he 
solved  automatically. 

Our  NCGB  software  is  a  long  complicated  C++,  Mathematica,  TeX  code  and  we  only 
support  it  on  Solaris  platforms.  Within  the  last  few  months,  the  C++  GB  engine  has 
been  heavily  modified  (modularized)  and  runs  under  Linux  and  Windows.  Thus  a  stand 
alone  C++  version  on  these  platforms  is  in  sight.  Testing  is  in  progress  as  is  linking  to 
Mathematica. 

Recently  we  (our  group  and  Bakshee  at  WRI)  found  a  way  to  give  Mathematica’s  control 
package  (by  Bakshee)  noncommuting  capability.  We  can  load  Bakshee’s  Mathematica  tool¬ 
box,  NCAlgebra,  and  our  new  Mathematica  file  file.m,  then  the  system  connection  laws 
portion  of  the  toolbox  works  with  noncommuting  A,B,C,D ’s.  Currently  we  are  extending 
the  file.m  to  make  those  parts  of  the  toolbox  (where  meaningful)  work  noncommutatively. 


This  should  lead  to  a  good  front  end  for  use  by  engineers  who  wish  to  conveniently  manip¬ 
ulate  many  noncommutative  linear  system  formulas  on  the  computer,  much  as  they  do  by 
hand. 

1.3  H°°  Design  and  Matrix  Inequalities 

1.3.1  H°°  Design 

This  concerns  a  basic  question  of  worst  case  frequency  domain  design  where  stability  of  the 
system  is  the  key  constraint.  This  is  the  H°°  optimization  problem  which  is  crucial  in  several 
branches  of  engineering. 

The  fundamental  if00  problem  of  control,  (with  0.  Merino  and  T.  Walker) 

First  we  state  the  core  mathematics  problem  graphically.  At  each  frequency  u  we  are  given 
a  set  Sw(c)  C  CN  The  objective  is  to  find  a  function  T  with  no  poles  in  the  R.H.P.  so  that 
each  T(juj)  belongs  to  Sw(c).  In  fact  there  is  a  simple  picture  to  think  of  in  connection  with 
a  design 


Figure  2 


Typically  there  is  a  nested  family  of  target  sets  Sw(c)  parameterized  by  a  performance  level 
c;  the  smaller  the  sets  the  better  the  performance.  For  the  optimal  c  a  solution  T  exists  but 
no  solution  exists  for  tighter  specs. 

The  Horowitz  templates  of  control  can  be  transformed  into  this  type  of  picture.  When  each 
Sw(c)  is  a  ” disk”  this  problem  is  solved  by  traditional  H°°  control.  This  graphical  problem 
can  be  formulated  analytically  in  terms  of  a  performance  function  T  as 

(OPT)  Given  a  positive  valued  function  T  on  R  x  CN  (which  is  a  performance 
measure), 

find  7*  >  0  and  /*  in  An  which  solve 

7*  =  inf  sup  r(u,/(jw))  . 

and  this  of  course  is  what  one  puts  in  a  computer.  Here  A jv  is  the  set  of  A-tuples  (/i, . . . ,  //v) 
of  RHP-stable  functions. 

Collaborators  and  I  have  a  very  broad  based  attack  on  the  problem  and  address  most  aspects 
of  it. 


From  qualitative  theory  to  numerical  algorithms  and  diagnostics. 

Work  under  the  contract  was 

(1)  to  link  the  theory  of  this  problem  with  semidefinite  programming.  (  In 
particular  to  improve  our  earlier  primal  dual  algorithms  using  some  ideas  from 
matrix  semidefinite  programming  ) 

(2)  to  begin  extension  of  our  OPT  theory  to  the  situation  where  the  performance 
T  contains  an  uncertain  parameter.  (This  gives  algorithms  and  diagnostics.) 

(3)  (with  G.  Balas)  we  have  worked  out  in  detail  the  synthesis  case  of  our 
diagnostics.  We  ran  a  few  cases  with  crude  codeing.  Results  were  intriguing, 
but  there  were  too  few  cases  to  know  if  they  were  at  all  representative.  Con¬ 
version  of  the  simplest  diagnostics  to  /x- tools  format  is  near  completion.  Then 
tests  will  be  easy  to  run.  Conventional  wisdom  about  optimization  is  that  the 
quality  of  software  and  the  users  enjoyment  are  considerably  aided  by  having 
good  diagnostics. 

Currently,  with  Harry  Dym  we  are  working  on  interior  point  methods  for  MIMO  H°° . 

1.3.2  Bi  H 00  control  or  Frequency  Response  Scheduling 

This  is  a  frequency  domain  type  of  gain  scheduling.  This  (rather  new  )  idea,  where  one 
schedules  on  frequency  response  functions,  leads  directly  to  canonical  problems  in  several 
complex  variables  for  which  a  fair  amount  of  theory  exists. 

An  example  is  H°°  control  for  a  plant  where  part  of  it  is  known  and  a  subsystem  8  is  not 
known,  that  is,  the  response  of  the  plant  at  “frequency”  s  is  P(s,8($)).  We  assume  that 
once  our  control  (closed  loop)  system  is  running,  we  can  identify  the  subsystem  8  on  line. 
Thus  the  problem  is  to  design  a  function  K  offline  that  uses  this  information  to  produce  a 
H 00  controller  via  the  formula  K(s,8(s)).  The  challenge  is  to  pick  I\  so  that  the  controller 
yields  a  closed  loop  system  with  H°°  gain  at  most  7,  no  matter  which  8  occurs.  While  this 
is  entirely  a  frequency  domain  problem,  it  has  the  flavor  of  gain  scheduling  and  one  might 
think  of  it  as  H°°  gain  scheduling.  Note  it  is  a  strict  analog  to  Linear  Parameter  Varying 
control  LPV  or  LFT  based  control,  approaches  currently  meeting  with  great  success. 

In  article  [Hprep]  ,  we  show  that  LMV  control  problems  are  equivalent  to  certain  problems  of 
interpolation  by  analytic  functions  in  several  complex  variables.  These  precisely  generalize 
the  classical  (one  complex  variable  )  interpolation  (AAK  -commutant  lifting)  problems  which 
lay  at  the  core  of  H°°  control.  These  problems  are  hard,  but  the  last  decade  has  seen 
substantial  success  on  them  in  the  operator  theory  community,  resulting  from  the  focus  of 
efforts  made  by  the  generation  of  mathematicians  who  followed  AAK-Nagy-Foias-Sarason. 

1.3.3  Matrix  Inequalities 

Also,  we  have  started  to  port  some  of  our  H°°  optimization  to  another  nonsmooth  problem, 
namely  matrix  inequalities.  One  such  result  concerns  coordinate  iteration  methods  of  opti¬ 
mization.  We  consider  optimization  of  the  largest  eigenvalue  of  a  smooth  selfadjoint  matrix 


valued  function  T(X,  Y)  of  two  vector  or  matrix  variables  X  and  Y.  A  typical  problem  one 
faces  in  control  design  are  matrix  versions  of  minimizing  in  Y  and  maximizing  in  X.  Also, 
minimizing  in  X  and  Y  is  an  important  problem.  We  shall  assume  that  T  is  concave  or 
convex  in  Y  and  separately  in  A’.  As  a  consequence,  assume  we  have  a  computer  package 
which  will  optimize  7  over  X  with  Y  frozen,  and  optimize  7  over  Y  with  X  frozen.  This 
is  realistic  when  the  joint  behavior  in  X  and  Y  is  bad,  since  to  use  existing  commercial 
software,  one  must  apply  it  to  each  coordinate  separately,  and  so  can  be  used  only  to  give  a 
“coordinate  iteration”  algorithm  . 

Result  On  “well  behaved  T”,  coordinate  iteration  always  gives  a  local  optimum  for  the  miny 
maxx  problem. 

Result  We  give  strong  evidence  (via  optimality  conditions)  that  coordinate  iteration  almost 
never  gives  a  local  solution  to  the  miny  minx  problem.  We  give  a  practical  test  to  tell  if  you 
have  hit  a  local  optimum. 


2  Technology  and  Interactions  1999 

2.1  Meetings 

Feb  97  School  of  Engineering  Gave  several  Seminars  UCSB  Kokotivic 
May  97  AFOSR  Contractors  Conf  Wright  Patterson  AFB  attended 


Sept  97 

Conference  "Operators  Systems  and  Linear  Algebra"  Kaiserslauten 
Univ  of  Kaiserslauten 

-  computer  algebra  for  engineering  systems 

Oct  97 

Pedagogical  Inst  of  Odessa 

-  computer  algebra  for  engineering  systems 

Dec  97  Conf  on  Decision  and  Control,  SanDiego 
Workshop  on  Nonlinear  H~\infty  control  -  opening  speaker 

-  Numerical  nonlinear  (jet  engine  application) 

-  Matrix  Optimization 


1998 

February  98  UCBerkeley  Math  Dept 
Analysis  seminar  -  Nonlinear  systems 


Operator  Theory  seminar  -  computer  algebra 


March  98  UCSanta  Barbara  -  Math  Dept  Colloquium  -  Nonlinear  Systems 

March  98  -  Regional  AMS  meeting  Kansas  Special  Session  in  Operator  Theory 
Non-linear  systems  and  operators 

AFOSR  Contractors  LA.  AFB  attended  May  98F0 


June  98 

International  Workshop  on  Operator  Theory  and  Appl. 


July  1998  MTNS 

-  Nonlinear  systems 

-  Computer  algebra  and  systems 

-  Optimization 


Special  Semester  on  Symbolic  Computation  at 
Math  Sci  Research  Inst ,  Berkeley 
"month  long  member"  -I  actually  stayed  2  weeks  over  3  visits, 
and  gave  one  talk  on  symbolic  calculation  for  engineering  problems. 


Dec  98  Conf  on  Decision  and  Control,  Orlando 

-  Nonlinear  Control 

-  Nonlinear  Control 

-  Numerical  nonlinear  (jet  engine  application) 

-  Matrix  Optimization 


1999 

June  1999  American  Conference  on  Control,  San  Diego 

—  Nonlinear  control 

—  Bi  Hinfinity  control 


AFOSR  Contractors  LA.  AFB 


attended 


May  98 


Dec  99  Conference  on  Decision  and  Control,  Pheonix 

—  Nonlinear  Control 

—  Nonlinear  Control 

—  Nonlinear  Control  (Robot) 

—  Computer  Algebra  in  Control 

—  $H“\infty$  Optimization 


3  Personnel  Supported 

SENIOR  PERSONNEL:  Helton,  Bill  Peter  Dower  PGR  Mark  Kennel  PGR  IV  Mikhail 

Shushick  VI  Trent  Walker  PGR  1 

Independent  Contractors:  Keller  Stankus  Merino 

GRADUATE  STUDENT  RESEARCHERS:  Dell  Kronewitter  GSR1  Juan  Camino  GSR  IV 
Mike  Hardt  GSR  IV  Joshua  Griffin  GSR  1  David  Glickenstien  GSR  1  Eric  Rowell  GSR  1 
Bob  Slobojidan  GSR  1  Jeremy  Martin  GSR  1  Maihong  Shi  GSR  3  Eric  Rowell  GSR  1 
UNDERGRAD  STUDENT  STAFF:  Marianna  Rauk  Campbell  Mike  Torre  Phillipe  Bergman 

3.1  Consultative  and  advisory 

The  main  US  Tokamak  is  in  San  Diego  at  General  Atomics  Corp.  and  my  former  student 
Mike  Walker  is  in  charge  of  installing  a  modern  control  system  on  it.  When  the  control  side 
of  things  is  moving  we  have  regular  discussions  about  what  to  do  next. 

3.2  Transitions 

We  maintain  and  expand  two  computer  programs  which  run  under  Mathematica  which  are 
publicly  available. 

NCAlgebra,  has  potential  applications  in  many  fields  (get  it  from  http://math.ucsd.edu/~ncalg) . 
This  is  the  established  package  which  gives  Mathematica  noncommuting  capability. 
OPTDesign,  our  classical  control  program  is  available  from  http :  //math  .ucsd .  edu/~helton. 
Our  book  [HMer98],  while  independent  of  the  software,  illustrates  its  use. 

Earlier  progress  reports  mentioned  discussions  on  antenna  design  with  NOSC  engineers  in 
San  Diego.  These  have  led  to  a  Navy  initiative  to  apply  H°°  techniques  developed  in  control 
problems  problems  to  antenna  matching  problems.  The  problem  is  challenging  because 
characteristics  (scattering  parameters)  of  an  antenna  (given  as  data  on  the  jui  axis  )  are 
complicated,  due  to  the  interaction  of  the  antenna  with  much  of  the  metal  on  the  ship.  We 
are  using  techniques  worked  out  by  Merino  and  I  under  previous  AFOSR  funding,  and  to 
some  extent  p  -  tools,  to  find  matching  circuits.  We  now  have  good  methods  to  determine 
if  existing  designs  are  far  from  the  theoretical  optimum.  Now  we  are  looking  at  low  order 
models.  Contact  at  Spawar:  Dr.  Dave  Schwartz  (619)  553  2021,  Jeff  Allen  (619)  553-6566. 


Ford  Motor  Co.  in  June  1999  gave  a  modest  gift  to  support  my  research.  This  is  the  second 
contribution  Ford  has  made  to  my  work.  The  possibility  of  applying  the  methodolgy  we 
developed  for  compressor  stall  control  to  engines  intrigues  them.  At  Ford:  Davor  Hrovat 
(313)  322-1492. 

Honors  and  Awards 

Guggenheim  Fellow  1985,  Outstanding  Paper  IEEE  Control  Society  1986 


4  Publications 

All  articles  which  appeared  in  1996. 

[HWee96]  J.  W.  Helton  and  F.  Weening:  “Some  systems  theorems  arising  from  the  Bieber- 
bach  conjecture,”  Journal  Nonlinear  and  Robust  Control,  vol  6,  (1996).  (misrecorded  as  94 
in  some  lists)  p  65-82 

[BH  1994,  1996]  J.  W.  Helton  and  J.  Ball:  “Viscosity  solutions  of  Hamilton-Jacobi  equations 
arising  in  nonlinear  i/oo-control.”  J.  Math  Systems  Estimation  and  Control  vol  6  num  1 
1996 

(JMSEC  published  an  announcement  1994,  put  the  paper  on  the  net,  in  1996  they  published 
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years  contributed  heavily  to  it.  The  patent  is  for  a  method  to  determine  the  ideal  limits  of 
an  antennas  power  consumption.  The  beginnings  of  the  work  were  my  papers  of  mine  in  the 
late  1970s.  Indeed  this  may  be  the  earliest  H°°  engineering  work  to  receive  a  patent. 


Part  II.  DISCLOSURE  OF  INVENTION  24  Nov  98 

1 .  GENERAL  PURPOSE : 

The  problem  of  broadband  impedance  matching  is 
maximizing  the  power  transfer  uniformly  over  a  frequency 
band  from  a  source  to  a  load. 

The  purpose  of  this  invention  is  to  compute  the  optimal  performance 
than  can  be  obtained  by  a  passive  or  lossless  matching  network. 


2 .  BACKGROUND : 

Engineers  utilize  a  variety  of  methods  for  constructing  broadband 
matching  networks.  Some  methods  may  attempt  to  maximize  power 
transfer  uniformly  over  a  frequency  band,  but  most  settle  for 
optimizing  at  a  selection  of  frequencies  within  the  band. 

These  methods  have  two  things  in  common: 

First,  the  ignorance  of  what  the  optimal  performance  can  be; 
Second,  the  lack  of  a  systematic  approach  to  optimize 
power  transfer  uniformly. 

This  invention  resolves  the  first  of  these  problems 
by  enabling  the  engineer  to  compute  the  optimal  performace 
for  a  matching  network  connecting  a  transmission  line  to 
a  load. 

The  application  that  motivated  this  effort  was  the 
need  to  maximize  power  radiated  by  HF  broadband  antennas 
whose  power  is  provided  from  a  50  Ohm  transmission  line. 

However,  this  approach  is  not  limited  to  any  frequency 
range  and  will  apply  to  microwave  and  electro-optic 
networks. 

Here  are  two  additional  benefits  of  knowledge  of  optimum 
performance  attainable  by  a  matching  circuit : 

First,  the  performance  optimum 

allows  an  engineer  a  benchmark  to  grade  the  performance 
of  proposed  matching  circuits. 

Second,  by  computing  the  optimum  performance  an  engineer 
can  determine  if  it  is  possible  to  satisfy  a  design  objective 
for  a  given  load. 


3.  DESCRIPTION  AND  OPERATION: 

Our  invention  is  an  algorithm  for  computing  the  performance  of  the 
optimum  matching  circuit  connecting  a  load  to  a  constant  impedance 
transmission  line  over  a  continuous  frequency  band  F. 

This  algorithm  is  implemented  in  a  computer  code  called  AIM  for 
Analytic  Impedance  Match. 

There  are  three  inputs  to  AIM: 

1.  qin  --  Samples  of  the  reflectance  of  the  load  normalized  to  the 
characteristic  impedance  of  the  line; 

2.  fin  --  The  frequencies  associated  with  the  normalized  reflectances; 

3.  N  --  A  parameter  specifying  the  number  of  frequencies  AIM  will 
compute  over . 

There  is  a  single  output  from  AIM: 

G  --  The  greatest  transducer  power  gain  obtainable  by  any  matching  circuit. 
There  are  two  sytem  parameters  computed  in  terms  of  G : 

1.  The  power  mismatch  (pm)  is  the  square  root  of  1-G; 

2.  The  voltage  standing  wave  ratio  VSWR  =  (1  +  pm)  /  (1  -  pm)  . 


.  NO 

I 


Thlf  algorithm  is  completely  developed  in  Attachment  1 . 

The  "language'  •  describing  this  algorithm  is  MATLAB  version  5.1  from 
The  Math  Works. 

The  main  program  is  the  "function1  1  AIM  with  the 
inputs  and  outputs  described  above. 

For  the  simplicity  of  this  implementation  N  is  assumed  to  be  a  power  of  2 . 

An  example  of  the  operation  of  AIM  follows: 

Here  qin  is  401  samples  from  a  data  set  labeled  AOC21B  and 
depicted  in  Figure  1. 

These  samples  correspond  to  frequencies  fin  equally  space  from 
2  to  8  MHz  inclusive. 

The  parameter  N  is  set  to  1024. 

The  output  of  AIM  is  G  =  0.705  as  the  largest  transducer 
power  gain  obtainable  by  a  matching  circuit . 

The  corresponding  power  mismatch  is  0.5431  and  VSWR  equal  to  3.378. 

4  .  ADVANTAGES  AND  NEW  FEATURES  : 

The  principle  advantage  of  this  invention  is  the  ability  to  compute  the 
optimal  performance  attainable  by  a  lossless  matching  network  connecting 
a  transmission  line  to  a  load  over  a  frequency  band. 

No  other  method  is  able  to  calculate  what  the  optimum  performance  is . 

Here  are  some  of  the  benefits  of  knowledge  of  optimum 
performance  attainable  by  a  matching  circuit : 

First,  the  ability  to  decide  whether  a  design  objective  for  a 
given  load  and  frequency  band  can  be  attained; 

Second,  a  benchmark  to  grade  the  performance 
of  proposed  matching  circuits . 


5 .  ALTERNATIVES : 

There  is  no  known  method  that  calculates  the  optimum  performance. 


6.  CONTRIBUTIONS  BY  INVENTORS: 

The  idea  for  this  method  was  developed  by  Professor  J.  William  Helton 
in  the  paper  "Broadbanding:  Gain  Equalization  Directly  from  Data'  ' 
published  in  IEEE  Transactions  on  Circuits  and  Systems  in  December 
1981.  The  implementation  was  done  by  David  F.  Schwartz  and  Jeffery  C. 
Allen. 

7.  EXECUTION  OF  DISCLOSURE: 


REPORT  DOCUMENTATION  PAGE 


Pubic  reporting  burden  for  mb  coiection  of  information  b  estimated  to  average  1  hour  per  response,  heading  me  time  for  revfa  A  rD  T  m  -rr\ 

gathering  and  maintalnkig  me  data  needed,  and  completing  and  reviewing  me  colecticn  of  information.  Send  comments  regarc  J\T  KL-^K-dL-  I K-UO- 
of  ^formation.  InckJdbig  suggestions  for  reducing  mb  burden  to  Washington  Headquarters  Service,  Directorate  for  trrformatfcn 

1215  Jefferson  Oavh  Highway,  Sute  1204,  Arlington,  VA  22202-4302.  and  to  me  Office  of  Management  and  3udget,  _  . 

Paperwork  Reduction  Project  (0704-0185)  Washetgton,  OC  20503.  f\  j 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS.  _  v-'V 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  DATE 

3/24/2000  _ 


4.  TITLE  AND  SUBTITLE 

Design  of  Robust  Controllers:  Frequency 
Domain  Methods  and  Their  Non-Linear 
Extensions 


4/1/97-12/31799 


1 5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


6.  AUTHOR(S) 


0.  William  Helton 


Award  #  F49620-97-0192 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUM3ER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONfTOR'S  ACRONYM(S) 

AFOSR 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  Mathematics 
Univ.  of  California,  San  Diego 
9500  Gilman  Drive 
La  Jolla,  CA  92093-0112 _ 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AFOSR 

110  Duncan  Avenue,  Ste  B115 
Bolling  AFB,  DC  20332-0001 


12.  DISTRIBUTION  AVAILABILITY  STATEMENT 

Approved  for  public  release, 
distribution  unlimited _ 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT  The  accomplishments  were  to  extend  H°°  control  in  new  directions: 

•  non-linear  control.  Theoretical  work  gave  algorithms  for  determining  in  real  time  key  properties  (the 
state)  of  a  system  one  tries  to  control  Purely  numerical  work  gave  a  "turnkey”  methodology  for  designing 
optimal  controllers  for  the  compressor  stall  problem.  Mathematically  related  was  numerical  solution  of  the 
extremely  hard  path  planning  problem,  biped  motion,  a  ” 5-link  robot”. 

•  .  computer  algebra  methods  appropriate  for  systems  which  are  aggregates  of  smaller  systems  (com¬ 
mercial  computer  algebra  packages  now  do  not  do  this). 

•  one  of  the  most  common  techniques  now  used  for  performance  optimization  (co-ordinate  descent 
applied  to  matrix  inequalities)  and  show  that  it  gives  the  wrong  answer  with  probablility  one.  Also  results 

_ _ give  a  test  to  see  how  far  off  (local  optimum)  this  is. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: _ 

a7  REPORT  I  b.  ABSTRACT  I  c.  THIS  PAGE 


17.  LIMITATION  OF 
ABSTRACT 


18.  NUMBER  19a. 
OF  PAGES 


NA*EOFRESppN^LEPERSOtN  . 

//.W,  \lw bh  7fn 

7a.  E^I^JIUMBE^p^ite  arMcode^^ 


EEIC  QUALITY  mS’BCTED  4 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI-Std  Z39-18 


